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Abstract 

A method for the semiclassical quantization of chaotic maps is proposed, which is 
based on harmonic inversion. The power of the technique is demonstrated for the 
baker's map as a prototype example of a chaotic map. 



The harmonic inversion method for signal processing [1,2] has proven to be a 
powerful tool for the semiclassical quantization of chaotic as well as integrable 
dynamical systems [3-5]. Starting from Gutzwiller's trace formula for chaotic 
systems, or the Berry- Tabor formula for integrable systems [6], the harmonic 
inversion method is able to circumvent the convergence problems of the pe- 
riodic orbit sums and to directly extract the semiclassical eigenvalues from a 
relatively small number of periodic orbits. The technique has successfully been 
applied to a large variety of Hamiltonian systems [4,5]. It has been shown that 
the method is universal in the sense that it does not depend on any special 
properties of the dynamical system. 

In this Letter we demonstrate that the range of application of the harmonic 
inversion method extends beyond Hamiltonian systems also to quantum maps. 
Starting from the analogue of Gutzwiller's trace formula for chaotic maps, we 
show that the semiclassical eigenvalues of chaotic maps can be determined 
by a procedure very similar to the one for flows. As an example system we 
consider the well known baker's map. For this map we can take advantage of 
the fact that the periodic orbit parameters can be determined analytically. 

We briefly review the basics of quantum maps that are relevant to what follows 
(for a detailed account of quantum maps see, e.g. Ref. [7]). We consider quan- 
tum maps, acting on a finite dimensional Hilbert space of dimension N, which 
possess a well-defined classical limit for iV — > oo. The quantum dynamics is 



Preprint submitted to Physics Letters A 



8 February 2008 



determined by the equation 



1p n+1 = Ulp n , (1) 



where U is a unitary matrix of dimension N, and ip n is the iV-dimensional 
discretized wave vector. The eigenvalues ut of U lie on the unit circle, u k = 
exp(—iip k ). The density of eigenphases ip k on the unit circle is given by 
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P qm (y?) = y + - Re £ Tr f/ n e"* , (2) 
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which can be rewritten as 



Aim(y?) = -- Im^ qm (v?) (3) 
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with the response function g qm given by 

gM=9M-^^ un ^ v ■ ( 4 ) 
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In analogy with the periodic orbit theory for flows, a semiclassical approxi- 
mation to the response function (4) can be obtained in terms of the periodic 
orbits of the corresponding classical system. [In the context of maps, "periodic 
orbit" means a sequence of fixed points periodic after n iterations; cyclic shifts 
of the same sequence correspond to the same periodic orbit.] In the semiclassi- 
cal approximation for maps, the traces of U n are related to the periodic orbits 
of topological length n, 



=■ Un , (5) 



Tr/7" ~ V — p i(Spo/n-Mpo-7r/2) _. 

~ P ^,|det(M po -l)|Vs e 

where the sum runs over all periodic orbits of topological length n including 
multiple traversals of shorter orbits. Here, S po is the action associated with 
the periodic orbit, /x po is its Maslov index, M po is the monodromy matrix of 
the orbit, and no is the topological length of the underlying primitive orbit 
(i.e., the length of the shortest subperiod). The value of the Planck constant is 
related to the dimension of the Hilbert space via H = l/(2nN). The accuracy 
of the semiclassical approximation can therefore be expected to improve with 
increasing N. 

The central idea for applying harmonic inversion to the periodic orbit quanti- 
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zation of maps now is to adjust the semiclassical response function 

<K*0 = <*>(¥>) + E^ e,B " ( 6 ) 



to the form of the exact quantum response function (4), expressed in terms of 
the eigenphases ipk and their multiplicities m^, 

<?qm(</?) = E ~~~~T ■ ( 7 ) 



It should be pointed out that for all maps by virtue of (4) the semiclassical 
amplitudes A n are independent of the phase ip. In analogy with the harmonic 
inversion procedures for Hamiltonian flows [3,4], we Fourier transform the 
oscillating part of the semiclassical response function (6) to obtain the semi- 
classical signal 

C(s) = J2A n 5(s-n). (8) 



The eigenphases ifk can now be determined by adjusting the semiclassical 
signal (8) to the form of the corresponding exact quantum signal (the Fourier 
transform of the exact response function (7)) 

Cqm(s) = -i£m fc e- u ^ (9) 
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by harmonic inversion. Note that compared to the corresponding procedure 
for flows, the topological length now plays the role of the scaled action, while 
the actions S po of the orbits are included in the amplitudes A n . Therefore, all 
orbits of the same topological length n contribute to the same signal point. 
While for flows the semiclassical signal takes on the simple form of a sum over 
5 functions only if we assume a certain scaling property [4,5], the form of the 
signal for maps is always the same, as the amplitudes A n are for all maps 
independent of ip. 

We will now apply the general procedure discussed above to the example of 
the baker's map, which has been used as a prototype example for studying 
the semiclassics of chaotic maps in many investigations in recent years [8-13]. 
The classical baker's map acts on points (p, q) of the unit square [0, 1] x [0, 1] 
according to 



q' = 2q mod 1 
p' = (p+[2q])/2, 



(10) 
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where [x] denotes the integer part of x. The periodic orbits of the map can be 
described by a complete binary symbolic code. Each orbit is characterized by 
a symbol string {ei, . . . , e„}, where q € {0, 1}. The action associated with a 
periodic orbit is given by [9] 



S„ = 



vv 



2 n - 1 



mod 1 



(12) 



with the integers v and v defined by 
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k=i 



n—k 



(13) 



Each orbit of length n has stability (i.e., largest eigenvalue of the monodromy 
matrix) 2 n . 

A quantized version of the baker's map which preserves the classical symme- 
tries was derived by Saraceno [8]. For even dimension N of the Hilbert space 
the quantum map is given by the unitary matrix 



U(N) = F^ 1 x 
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'N/2 

F, 







(14) 
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with 



{F N )nm = -J= exp 



N 



- 2l d{n-\){m-\) 



n,m= 1,...,N. (15) 



We will use the the exact quantum eigenvalues for comparison with the semi- 
classical ones obtained by harmonic inversion of the semiclassical signal (8). 

For the baker's map, the semiclassical amplitudes A n in the signal (8) as 
defined in Eq. (5) read [11] 

A-iE^TT^M^TT (16) 

po(n) 



with v and v given by Eq. (13). We have performed calculations for dimensions 
iV = 6 and N — 12. The case N = 6 has also been examined in Ref. [11], where 
the semiclassical eigenvalues were determined by a resummation of the Selberg 
zeta function of the map. For the construction of the semiclassical signal, we 
calculated all orbits up to symbol length n = 20 for dimension N = 6 and up 
to length n = 38 for iV = 12. Figure 1 shows the results for the eigenvalues 
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Uk = exp(— i(pk) obtained by harmonic inversion of the semiclassical signal (□), 
compared with the exact quantum results (*) obtained by diagonalization of 
the quantum matrix U (see Eq. (14)). For comparison, in the case N = 6, we 
have also plotted the semiclassical results from Ref. [11] (+). 

Our semiclassical results are in good agreement with the exact quantum eigen- 
values. However, as was also pointed out in Ref. [11], in the case of the baker's 
map the semiclassical error is relatively large, which is probably due to the 
discontinuities inherent in this map. In particular, a few of the semiclassical 
eigenvalues are located away from the unit circle at distances on the order 
of KT 1 . Although the accuracy of the semiclassical approximation should im- 
prove in general with increasing N, the doubling of the dimension of the 
Hilbert space from iV = 6toiV = 12on average does not yet visibly lead the 
semiclassical eigenvalues closer to the unit circle. 

On the other hand, in the case N = 6, the semiclassical eigenvalues obtained by 
harmonic inversion are in excellent agreement with those from Ref. [11]. This 
implies that the deviations from the exact quantum eigenvalues are indeed 
solely due to the semiclassical error, and do not indicate any inaccuracies of 
the individual methods applied. We add that the semiclassical approximation 
for the baker's map may even be improved by introducing non-semiclassical 
correction factors to the amplitudes A n [11,12]. 

As for flows, the harmonic inversion method for maps can also be used, vice 
versa, to analyze the quantum spectrum in terms of the amplitudes A n [14]. 
This is achieved by adjusting the quantum response function (7) to its semi- 
classical approximation (6). However, since for maps the amplitudes A n con- 
tain contributions from all periodic orbits with topological length n, it is not 
possible to extract information about single periodic orbits from the quantum 
spectrum. Moreover, since the traces TrU n do not depend on the quantity 
ip which is to be quantized, the exact expression (4) for maps already fulfills 
the ansatz of the harmonic inversion procedure as a whole. This is in contrast 
to scaling Hamiltonian systems, where the higher order h corrections to the 
semiclassical approximation depend on the scaling parameter (e.g., for billiard 
systems, the wave number) and the harmonic inversion of the exact quantum 
spectrum yields only the zeroth order H contributions to the response function, 
with the higher orders acting as a kind of noise [4,5]. As a consequence, the 
analysis of the quantum spectrum of maps will simply yield the exact quan- 
tum values for the traces Tr U n rather than their semiclassical approximations 
A n (we have verified this for the baker's map), and no information about the 
semiclassics or single periodic orbits can be obtained. 

In conclusion, we have presented a method for the periodic orbit quantization 
of chaotic maps, which makes use of harmonic inversion. The procedure works 
similar to the one for flows, and can in the same way be applied to all chaotic 
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maps, independent of any special properties of the respective system. We 
have demonstrated the power of the method by successfully applying it to 
the baker's map. Our results reproduce the exact quantum eigenvalues of the 
baker's map to within the error of the semiclassical approximation, and are in 
excellent agreement with those obtained by other semiclassical methods. 

We gratefully acknowledge support of this work by the Deutsche Forschungs- 
gemeinschaft. 
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Fig. 1. Exact quantum (*) and semiclassical (□) eigenvalues Uk = exp(— \<pk) of 
the baker's map obtained by harmonic inversion for dimensions (a) iV = 6 and (b) 
N = 12. In the case N = 6 the semiclassical results from Ref. [11] are also plotted 



